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ABSTRACT 

When our Sun was young it rotated much more rapidly than now. Observations of young, rapidly rotating 
stars indicate that many possess substantial magnetic activity and strong axisymmetric magnetic fields. We 
conduct simulations of dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic 
(ASH) code to explore the complex coupling between rotation, convection and magnetism. Here we study 
dynamo action realized in the bulk of the convection zone for a system rotating at three times the current solar 
rotation rate. We find that substantial organized global-scale magnetic fields are achieved by dynamo action in 
this system. Striking wreaths of magnetism are built in the midst of the convection zone, coexisting with the 
turbulent convection. This is a surprise, for it has been widely believed that such magnetic structures should be 
disrupted by magnetic buoyancy or turbulent pumping. Thus, many solar dynamo theories have suggested that 
a tachocline of penetration and shear at the base of the convection zone is a crucial ingredient for organized 
dynamo action, whereas these simulations do not include such tachoclines. We examine how these persistent 
magnetic wreaths are maintained by dynamo processes and explore whether a classical mean-field a-effect 
explains the regeneration of poloidal field. We find that the global-scale toroidal magnetic fields are maintained 
by an il-effect arising from the differential rotation, while the global-scale poloidal fields arise from turbulent 
correlations between the convective flows and magnetic fields. These correlations are not well represented by 
an a-effect that is based on the kinetic and magnetic helicities. 

Subject headings: convection - MHD - stars interiors - stars:rotation - stars: magnetic fields - Sun:interior 



1 . STELLAR MAGNETISM AND ROTATION 

Most stars are born rotating quite rapidly. They can ar- 
rive on the main sequence with rot ational velocities as high 
as 200 km s" 1 dBouvier et all 119971) . Stars with convec- 
tion zones at their surfaces, like the Sun, slowly spin down 
as they shed angul ar momentum through their magnetized 
stellar winds (e.g.. IWeber & DavisHl967t ISkumanichlll972l; 
iMacGregor & Brennerlll99lTL The time needed for signifi- 
cant s pindown appears to be a stron g function of stellar mass 
(e.g., [Barnes 20 031 : IWest et al.ll2004ft : solar-mass stars slow 
less rapidly than somewhat less massive G and K-type stars, 
but still appear to lose much of their angular momentum by 
the time they are as old as the Sun. Present day observations 
of the solar wind likewise indicate that the current angular 
momentum flux from the Sun is a few times 10 30 dyn cm (e.g., 
iPizzo et aDl 19831) . suggesting a time scale for substantial an- 
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gular momentum loss of a few billion years. Thus the Sun 
likely rotated significantly more rapidly in its youth than it 
does today. 

1.1. Rotation-Activity Relations 

Rotation appears to be inextricably linked to stellar mag- 
netic activity. Observations indicate that in stars with ex- 
tensive convective envelopes, chromospheric and coronal ac- 
tivity - which partly trace magnetic heating of stellar atmo- 
spheres - first rise with increasing rotation rate, then even- 
tually level off at a constant value for rotati on rates above a 
mass-dependent threshold velocity (e . g.. iNoyes et al] 119841; 
[ Patten & Simon] 119961: iDelfosse et at] il998t IPizzolato et all 
2003). Activit y may even decline somewhat in the most rapid 
rotators (e.g., James et al] 120001) . Similar correspondence 
is observed between rotation rate and estimates of the un- 
signe d surface magnetic flux (Saar 19961 l200lt iReiners et alj 
2009). This "rotation-activity" relationship is tightened when 
stellar rotation is given in terms of the Rossby number 
Ro ~ PJt c , with P the rotation period and r c an estimate 
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of the convective overturning time (e.g.. iNoves et al.lll984l) . 
Expressed in this fashion, a common rotation-activity cor- 
relation appears to span spectral types ranging from late 
F to late M (e.g.. iMohanty & Basril l200l iPizzolato et alJ 
12001 iReiners & Basril 12007b . Magnetic fields can likewise 
feed back upon stellar rotation by modifying the rate at 
which angular momen t um is lost through a st ellar wind (e.g., 
IWeber & Davisl 119671; iMatt & Pudritzl 120081) . Analyses of 
stellar spindown as a function of age and mass have thus pro- 
vided further constraints on stellar magnetism and its connec- 
tions to rotation. There are also indications that the period 
of the activity cycle itself may depend o n the stellar rota- 
tion rate (e.g.. ISaar & Bra ndenburg 1999). Recent observa- 
tions of solar-type stars may indicate that even the topology 
of the global-scale fields changes with rotation rate, with the 
rapid rotators having substan tial global-scale toroidal mag- 
netic fields at their surfaces (Petit et al. 2008). The overall 
picture that emerges from these observations is that rapid ro- 
tation, as realized in the younger Sun and in a host of other 
stars, can aid in the generation of strong magnetic fields and 
that young stars tend to be rapidly rotating and magnetically 
active, whereas older ones are slow er and less active (e.g., 
lBarnesll2003l: IWest et al.ll200l 120081) . 

A full theoretical understanding of the rotation-activity re- 
lationship, and likewise of stellar spindown, has remained 
elusive. Some aspects of these phenomena probably de- 
pend upon the details of magnetic flux emergence, chro- 
mospheric and coronal heating, and mass loss mechanisms 
- but the basic existence of a rotation-activity relationship 
is widely thought to reflect some underlyin g rotational de- 
pendence of the dynamo process itself (e.g., iKnobloch et alJ 
fl98lt INoves el^ll9MlBaliunas et al.lfl996h . 

1.2. Elements of Global Dynamo Action 

In stars like the Sun, the global-scale dynamo is gener- 
ally thought to be seated in the tachocline, an interface of 
shear between the differentially rotating convection zone and 
the radiative interior which is in solid body rotation (e.g., 
Parkerill993l ICharbonneau & MacGregoril 19971; Ibssendriiveri 
20031) . Helioseismology revealed the internal rotation profile 
of the Sun and the pres ence of this important shear layer (e.g., 
Thom pson et al.l 120031) . The stably stratified tachocline may 
provide a region for storing and amplifying coherent tubes 
of magnetic field which may eventually rise to the surface of 
the Sun as sunspots. Others have suggested that the latitu- 
dinal and radial gradients of angular velocity in the bulk of 
the convec tion zone may be sufficient for global dynamo ac- 
tion (e.g., [D ikpati & Charbonneau 1999; Brandenburgi l2005l; 
iGuerrero & de Gouveia Dal Pinol l2007). However, it has gen- 
erally been believed that magnetic buoyancy instabilities may 
prevent fields from being s trongly ampl ified within the bulk of 
the convection zone itself (Par ker! 19751) . In the now prevalent 
"interface dynamo" model, solar magnetic fields are partly 
generated in the convection zone by helical convection, then 
transported downward into the tachocline where they are or- 
ganized and amplified by the shear. Ultimately the fields may 
become unstable and rise to the surface. 

Although the rotational dependence of this process is not 
well understood, some guidance may come from mean-field 
dynamo theory. In such theories, the solar dynamo is often 
referred to as an "a - 17" dynamo, with the a-effect ch aracter- 
izing t he twisting of fie lds by helical convection (e.g., Moffatt 
119781: Steenbec k et alj[l966l) . and the Sl-effect representing 
the shearing of poloidal fields by differential rotation to form 



toroidal fields. Both of these effects are, in mean-field theory, 
sensitive to rotation: the a-effect because it is proportional to 
the kinetic helicity of the convective flows, which sense the 
overall rotation rate, and the £!-effect because more rapidly 
rotating stars are generally expected to have stronger differ- 
ential rotation. But the detailed nature of these effects in the 
solar dynamo and the appropriate scaling with rotation has 
been very difficult to elucidate. 

Simulations of the global-scale solar dynamo have gener- 
ally affirmed the view that the tachocline may play a central 
role in building the globally-ordered magnetism in the Sun. 
Recent three-dimensional (3D) simulations of solar convec- 
tion without a tachocline at the base of the convection zone 
achieved dynamo action and produced magnetic fields which 
were strongly domin ated by fluctuatin g components with lit- 
tle global-scale order (Br un et al.l2004l) . When a tachocline of 
penetration and shear was included, remarkable global-scale 
magnetic structures were realized in the tachocline region, 
while the convection zone re mained dominated by fluctuating 
fields (Browning et al. 2006). These simulations are making 
good progress toward clarifying the elements at work in the 
operation of the solar global-scale dynamo, but for other stars 
many questions remain. In particular, observati ons of large- 
scale magnetism in fully convective M-stars (Don atiet alJ 
2006), along with the persistence of a rotation-activity cor- 
relation in such low-mass stars, hint that perhaps tachoclines 
may not be essential for the generation of global-scale mag- 
netic fields. This view is partly borne out by simulations of M - 
dwarfs under strong rotational constraints dBrown ing 2008), 
where strong longitudinal mean fields were realized despite 
the lack of either substantial differential rotation or a stable in- 
terior and thus no classical tachocline. Major puzzles remain 
in the quest to understand stellar magnetism and its scaling 
with stellar rotation. 

1.3. Convection and Dynamos in Rapidly Rotating Systems 

We began our study of rapidly rotating suns by carrying 
out a suite of 3D hydrodynamic simulations in full spherical 
shells that explored the coupling o f rotation and co nvection 
in these younger solar-type stars dBrown et alJl2008h . Those 
simulations studied the influence of rotation on the patterns of 
convection and the nature of global-scale flows in such stars. 
The shearing flows of differential rotation generally grow in 
amplitude with more rapid rotation, possessing rapid equators 
and slower poles, while the meridional circulations weaken 
and break up into multiple cells in radius and latitude. More 
rapid rotation can also substantially modify the patterns of 
convection in a surprising fashion. With more rapid rotation, 
localized states begin to appear in which the convection at 
low latitudes is modulated in its strength with longitude. At 
the highest rotation rates, the convection can become confined 
to active nests which propagate at distinct rates and persist for 
long epochs. 

Motivated by these discoveries, we turn here to explorations 
of the possible dynamo action achieved in a solar-type star ro- 
tating at three times the current solar rate. These 3D magneto- 
hydrodynamic (MHD) simulations span the convection zone 
alone, as the nature of tachoclines in more rapidly rotating 
suns is at present unclear. We find that a variety of dynamos 
can be excited, including steady and oscillating states, and 
that dynamo action is substantially easier to achieve at these 
faster rotation rates than in the solar simulations. Magnetism 
leads to strong feedbacks on the flows, particularly modify- 
ing the differential rotation and its scaling with the overall 



Persistent Wreaths of Magnetism 



3 



rotation rate fi(j. The magnetic fields which form in these 
dynamos have prominent global-scale organization within the 
convect ion zone, in contrast to previous solar dynamo simu- 
lations dBrun et al.ll2004tlBrowning et al.ll2006l) . 

Quite strikingly, we find that coherent global magnetic 
structures arise naturally in the midst of the turbulent convec- 
tion zone. These wreath-like structures are regions of strong 
longitudinal field organized loosely into tubes, with fields 
wandering in and out of the surrounding convection. These 
wreaths of magnetism differ substantially from the idealized 
flux tubes supposed in many dynamo theories, though they 
may be related to coherent structures achie ved in local sim- 
ulations of dynamo action in shear flows dCline et al.ll2003t 
iVasil & Brummellll2008ll2009h . 

Here we explore the nature of persistent magnetic wreaths 
realized in a global simulation rotating at three times the so- 
lar rotation rate, and discuss how they are maintained amidst 
turbulent convection. In many of our other rapidly rotat- 
ing suns, the dynamos become time dependent and undergo 
semi -regular changes of global-scale polarity. Those dynamos 
will be explored in an upcoming paper. We additionally find 
that magnetic wreaths survive in the presence of a model 
tachocline, and those simulations will be reported on sepa- 
rately. 

We outline in ^2] the 3D MHD anelastic spherical shell 
model and the parameter space explored by these simulations. 
We then examine in §^3]and|4]the structure of magnetic fields 
found in our rapidly rotating dynamo at three times the so- 
lar rate, which builds persistent global-scale ordered fields in 
the form of wreaths in the midst of its convection zone. In 
$5] we examine how such global-scale fields are created and 
maintained by dynamo processes. In $6] we explore whether 
a classical mean-field a-effect reproduces our observed pro- 
duction of poloidal field. We reflect on our findings in $7] 

2. GLOBAL MODELLING APPROACH 

To study the coupling between rotation, magnetism and 
the large-scale flows achieved in stellar convection zones, we 
must employ a global model which simultaneously captures 
the spherical shell geometry and admits the possibility of 
zonal jets and large eddy vortices, and of convective plumes 
that may span the depth of the convection zone. The solar 
convection zone is intensely turbulent and microscopic values 
of viscosity and magnetic and thermal diffusivities in the Sun 
are estimated to be very small. Numerical simulations cannot 
hope to resolve all scales of motion present in real stellar con- 
vection and must instead strike a compromise between resolv- 
ing dynamics on small scales and capturing the connectivity 
and geometry of the global scales. Here we focus on the latter 
by studying a full spherical shell of convection. 

2.1. Anelastic MHD Formulation 

Our tool for exploring MHD stellar convection is the 
anelastic s pherical harmonic ( ASH) code, which is described 
in detail in lClune et al.l (11999b. The imp lementation of mag- 
netism is discussed in lBrun et al.l d2004l) . ASH solves the 3D 
MHD anelastic equations of motion in a rotating spherical 
shell using the pseudo-spectral method and runs efficiently on 
massively parallel architectures. We use the anelastic approx- 
imation to capture the effects of density stratification with- 
out having to resolve sound waves which have short periods 
(about 5 minutes) relative to the dynamical time scales of the 
global scale convection (weeks to months) or possible cycles 
of stellar activity (years to decades). This criteria effectively 



filters out the fast magneto-acoustic modes while retaining the 
slow modes and Alfven waves. Under the anelastic approxi- 
mation the thermodynamic fluctuating variables are linearized 
about their spherically symmetric and evolving mean state, 
with radially varying density p, pressure P, temperature T and 
specific entropy S. The fluctuations about this mean state are 
denoted as p, P, T and S. In the reference frame of the star, 
rotating at average rotation rate £Iq, the resulting MHD equa- 
tions are: 

V-(,ot) = 0, 



dv 
dt 



VB = Q, 



+ (v- V)u + 2fi X v 
1 



= -V(P+P) 



+(P+P)9+ — (VxB)xB-V-V, 

4-TT 



dB 



Vx(»x 6)-Vx (i]V X B), 
~dS 



(1) 
(2) 

(3) 
(4) 



pT 



Of 



+ v- V(S+S) 



V- [K r pc p V(T + T) + K Q pTVS+KpTVS] (5) 



j + 2pv 



1 



where v = (v^ve,^) is the local velocity in the stellar ref- 
erence frame, B = {B ri Be 1 B ( f,) is the magnetic field, j is the 
vector current density, g is the gravitational acceleration, c p is 
the specific heat at constant pressure, K r is the radiative diffu- 
sivity and T> is the viscous stress tensor, given by 



T>ij = —2pv 



(6) 



where e,j is the strain rate tensor. Here v, k and 77 are the dif- 
fusivities for vorticity, entropy and magnetic field. We assume 
an ideal gas law 

P = Tlpf, (7) 

where 1Z is the gas constant, and close this set of equations 
using the linearized relations for the thermodynamic fluctua- 
tions of 

p _P T _ P S 
p~P~T~^P~~c~ p 



(8) 



The mean state thermodynamic variables that vary with radius 
are evolved with the fluctuations, thus allowing the convection 
to modify the entropy gradients which drive it. 

The mass flux and the magnetic field are represented with a 
toroidal-poloidal decomposition as 

pv = Vx Vx(Wr) + Vx(Zr), (9) 
B=VxVx(/?r) + Vx((r), (10) 

with streamfunctions W and Z and magnetic potentials (3 
and This approach ensures that both quantities remain 
divergence-free to machine precision throughout the simula- 
tion. The velocity, magnetic and thermodynamic variables 
are all expanded in spherical harmonics for their horizon- 
tal structure and in Chebyshev polynomials for their radial 
structure. The solution is time evolved with a second-order 
Adams-Bashforth/Crank-Nicolson technique. 
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TABLE 1 

Parameters for Primary Simulations 



Case N r ,N e ,N^ Ra Ta Re Re' Rm Rm' Ro Roc v t] fl /fl Q 

D3 96x256x512 3.22xl0 5 1.22xl0 7 173 105 86 52 0.378 0.311 1.32 2.64 3 
H3 96x256x512 4.10xl0 5 1.22xl0 7 335 105 — — 0.427 0.353 1.32 — 3 



NOTE. — Dynamo simulation at three times the solar rotation rate is case D3, and the hydrodynamic (non-magnetic) 
companion is H3. Both simulations have inner radius n, 0l = 5.0 X 10 10 cm and outer radius of r t0 p = 6.72 X 10 10 cm, with 
L = (r to p — r^ot) = 1 .72 X 10 10 cm the thickness of the spherical shell. Evaluated at mid-depth are the Rayleigh number 
Ra = (-dp/dS)(dS/dr)gL 4 /pvK, the Taylor number Ta = 4f2 2 L 4 /^ 2 , the rms Reynolds number Re = v^^Ljv and fluctu- 
ating Reynolds number Re' = v' ms L/v, the magnetic Reynolds number Rm = v rms L/»7 and fluctuating magnetic Reynolds 
number Rm' = v' rms L/rj, the Rossby number Ro = oj/2f7o . and the convective Rossby number Roc = (Ra/TaPr)'/ 2 . Here 
the fluctuating velocity v' has the axisymmetric component removed: v' = v — (v), with angle brackets denoting an average 
in longitude. For both simulations, the Prandtl number Pr = vjn is 0.25 and in the dynamo simulation the magnetic Prandtl 
number Pm = i//rj is 0.5. The viscous and magnetic diffusivity, v and rj, are quoted at mid-depth (in units of 10 12 cm 2 s _1 ). 
The rotation rate f2o of each reference frame is in multiples of the solar rate Q© = 2.6 X 10 -6 rad s _1 or 414 nHz. The viscous 
time scale at mid-depth t v = L 2 /u is about 2600 days for case D3 and the resistive time scale is about 1300 days, while the 
rotation period is 9.3 days. 



ASH is a large-eddy simulation (LES) code, with subgrid- 
scale (SGS) treatments for scales of motion which fall be- 
low the spatial resolution in our simulations. We treat these 
scales with effective eddy diffusivities, v, k and 77, which 
represent the transport of momentum, entropy and magnetic 
field by unresolved motions in the simulations. These sim- 
ulations are based on the hydrodynamic studies reported in 
iBrown et al.l (|2008), and as there v, n and r\ are taken for 
simplicity as functions of radius alone and proportional to 
p~ l l\ This adopted SGS variation, as in iBrun et all (120041) 
and Browning et al. (2006), yields lower diffusivities near the 
bottom of the layer and thus higher Reynolds numbers. Act- 
ing on the mean entropy gradient is the eddy thermal diffusion 
kq which is treated separately and occupies a narrow region 
in the upper convection zone. Its purpose is to transport en- 
tropy through the outer surface where radial convective mo- 
tions vanish. 

The boundary conditions imposed at the top and bottom of 
the convective unstable shell are: 

1 . Impenetrable top and bottom: v r = , 

2. Stress-free top and bottom: 

(d/dr)(vg/r) = (d/dr)( V4> /r) = 0, 

3. Constant entropy gradient at top and bottom: 

d(S +S)/dr= const, (11) 

4. Match to external potential field at top: 

B = V<S> and V 2 $ = 0|,._, , 

1 ' — ' lop 

5. Perfect conductor at bottom: 

B r = (d/dr)(rB e ) = {d/dr){rB^ = . 

2.2. Posing the Dynamo Problem 

Our simulations are a simplified picture of the vastly tur- 
bulent stellar convection zones present in G-type stars. We 
take solar values for the input entropy flux, mass and radius, 
and explore simulations of a star rotating at three times the 
current solar rotation rate. We focus here on the bulk of the 
convection zone, with our computational domain extending 



from Q.IIRq to O.97/? , thus spanning 172 Mm in radius. 
The total density contrast across the shell is about 25. The 
reference or mean state of our thermody namic variables i s de- 
rived from a ID solar structure model (Brun et al. 2003) and 
is continuously updated with the spherically symmetric com- 
ponents of the thermodynamic fluctuations as the simulations 
proceed. The referen ce state in all of thes e simulations is sim- 
ilar to that shown in [Brown et al.l (120081) . We avoid regions 
near the stellar surface where hydrogen ionization and radia- 
tive losses drive intense convection (like granulation) on very 
small scales that we cannot resolve, and thus position the up- 
per boundary slightly below this region. Our lower boundary 
is positioned near the base of the convection zone, thus omit- 
ting the stably stratified radiative interior and the shear layer at 
the base of the convection zone known as the tachocline. The 
fundamental characteristics of our simulations and parameter 
definitions are summarized in Table Q] 

The dynamo simulation was initiated from a mature hydro- 
dynamic progenitor which had been evolved for more than 
5000 days and was well equilibrated. The progeni t or cas e H3 
is very similar to case G3 reported in Brown et al. (20Q8|), but 
here we chose a functional form for the SGS entropy diffusion 
kq that is more confined to the upper 10% of the convection 
zone; the unresolved flux here does not vary as much with 
rotation rate. The effects of this change are subtle, resulting 
primarily in slightly stronger latitudinal gradients of differen- 
tial rotation and temperature in the uppermost regions of the 
shell. The patterns of convection are very similar to those 
found in case G3, though here they are slightly more complex 
near the top of the shell, and the Reynolds number remains 
high throughout the convection zone. Case H3 possesses in- 
tricate convective patterns and a solar-like differential rotation 
profile, with fast zonal flow at the equator and slower flows at 
the poles. 

To initiate our dynamo case, a small seed dipole magnetic 
field was introduced and evolved via the induction equation. 
The energy in the magnetic fields is initially many orders of 
magnitude smaller than the energy contained in the convective 
motions, but these fields are amplified by shear and grow to 
become comparable in energy to the convective motions. 

Stellar dynamo simulations are computationally intensive, 
requiring both high resolutions to correctly represent the ve- 
locity fields and long time evolution to capture the equili- 
brated dynamo behavior, which may include cyclic variations 
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FIG. 1. — Convective structures and mean flows in cases D3 and H3. (a) Radial velocity v r in dynamo case D3, shown in global Mollweide projection at 
O.95R0, with upflows light and downflows dark. Poles are at top and bottom and the equator is the thick dashed line. The stellar surface at Rq is indicated by 
the thin surrounding line, (b) Profiles of mean angular velocity Q(r, 8), accompanied in (c) by radial cuts of Q at selected latitudes. A strong differential rotation 
is established by the convection, (d) Profiles of meridional circulation, with sense of circulation indicated by color (red counter-clockwise, blue clockwise) and 
streamlines of mass flux overlaid, (e-h) Companion presentation of fields for hydrodynamic progenitor case H3. The patterns of radial velocity are very similar 
in both cases. The differential rotation is much stronger in the hydrodynamic case and the meridional circulations there are somewhat weaker, though their 
structure remains similar. 



on time scales of several years. The strong magnetic fields can 
produce rapidly moving Alfven waves which seriously restrict 
the Courant-Friedrichs-Lewy (CFL) timestep limits in the up- 
per portions of the convection zone. Case D3, rotating three 
times faster than the current Sun, has been evolved for over 
7000 days (or over 2 million timesteps). We plan to report 
on a variety of other dynamo cases, some at higher turbulence 
levels and rotation rates, in subsequent papers. 

This dynamo simulation was conducted at magnetic Prandtl 
number Pm = v/rj = 0.5, a value significantly lower than 
emplo yed in ou r prev ious solar simulations. In partic- 
ular. IBrun et alj (120041) explored Pm = 2,2.5 and 4, and 
iBrowning et alj (120061) studied Pm = 8. The high magnetic 
Prandtl numbers were required in the solar simulations to 
reach sufficiently high magnetic Reynolds numb ers to drive 
sustain ed dynamo action. In the simulations of IBrun et alj 
(2004) only the simulations with Pm > 2.5 and Rm' > 300 
achieved sustained dynamo action, where Rm' is the fluctu- 
ating magnetic Reynolds number. We are here able to use 
a lower magnetic Prandtl number for three reasons. Firstly, 
more rapid rotation tends to stabilize convection and lower 
values of v and r\ are required to drive the convection. Once 
convective motions begin, they become quite vigorous and 
the fluctuating velocities saturate at values comparable to our 
solar cases. Thus the Reynolds numbers achieved are fairly 
large and we can achieve modestly high magnetic Reynolds 
numbers even at low Pm. Secondly, the differential rotation 
becomes substantially stronger with both more rapid rotation 
flo and with lower diffusivities v and r\. This global-scale flow 
is an important ingredient and reservoir of energy for these 
dynamos, and the increase in its amplitude means that low 



Pm dynamos can still achieve large magnetic Reynolds num- 
bers based on this zonal flow. Lastly, the critical magnetic 
Reynolds number for dynamo action likely decrea ses with in- 
creasing kinetic helicity (e.g. JLeorat et al.lll981|). and helic- 
ity ge nerally increases with rotation rate (e.g.. iKapvla et all 
2009). Indeed there are even suggestions that the presence 
of a mean shearin g flow may lower the c ritical magnetic 
Reynolds number (Hughes & Proctor 2009), and the strong 
differential rotation present in these rapidly rotating suns may 
serve to lower this threshold for dynamo action. We find that 
the rapidly rotating flows considered here achieve dynamo ac- 
tion a t somewhat lower Rm than the models of IBrun et alj 
(2004), which rotated at the solar rate. 

3. DYNAMOS WITH PERSISTENT MAGNETIC WREATHS 

We here explore case D3 which yields fairly persistent 
wreaths of magnetism in its two hemispheres, though these 
do wax and wane somewhat in strength once established. Ex- 
amining the properties of this dynamo solution should help to 
provide a perspective for the greater variations realized in our 
time-dependent dynamos which will be discussed in a follow- 
ing paper. 

3.1. Patterns of Convection 

The complex and evolving convective structures in our dy- 
namo cases are substantially similar to the patterns of con- 
vection found in our hydrodynamic simulations. Our dynamo 
solution rotating at three times the solar rate, case D3, is pre- 
sented in Figure [T] along with its hydrodynamic progenitor, 
case H3. The radial velocities shown near the top of the sim- 
ulated domain (Figs. QJ[,<?) have broad upflows and narrow 
downflows as a consequence of the compressible motions. 
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FIG. 2. — Magnetic wreaths and convective flows sampled at the same instant in case D3. (a) Longitudinal magnetic field B^ near the top of the shell (O.95S0) 
and (b) at mid-depth (O.85/?0). Strong flux structures with opposite polarity lie above and below the equator and span the convection zone. (c,d) Weaker radial 
magnetic field B, permeates and encircles each wreath, (e, /) Strong convective upflows and downflows shown by V r pass through and around the wreaths. The 
regions of strong magnetism tend to disrupt the convective flows while the strongest downflows serve to pump the wreaths to greater depths. 

in part to the thermal structure of th e solar tachocli ne, as first 
found in the mean-field models of Rempel (2005) and then 
in 3D simulations of global-scale convection by Miesc h et al.1 
(2006). In those computations, it was realized that introduc- 
ing a weak latitudinal gradient of entropy at the base of the 
convection zone, consistent with a thermal wind balance in 
a tachocline of shear, can serve to tilt the ft contours to- 
ward a more radial alignment without significantly chang- 
ing either the overall £1 contrast with latitude or the convec- 
tive patterns. iBallot et ail (120071) explored the consequences 
of such a boundary condition in one of their simulations of 
young, rapidly rotating suns with deep convection zones and 
found that the effect s on the differentia l rotation were simi- 
lar to those found in Miesch et al. (2006). We expect similar 
behavior here, but at present observations of rapidly rotating 
stars only measure differential rotation at the surface and do 
not offer constraints on either the existence of tachoclines in 
young suns or the nature of their internal differential rotation 
profiles. As such, we have neglected the possible tachoclines 
of penetration and shear entirely in these models and instead 
adopt the simplification of imposing a constant radial entropy 
gradient at the bottom of the convection zone. 

The differential rotation achieved is stronger in our hydro- 
dynamic case H3 than in our dynamo case D3. This can be 



Near the equator the convection is aligned largely in the north- 
south direction, and these broad fronts sweep through the do- 
main in a prograde fashion. The strongest downflows pene- 
trate to the bottom of the convection zone; the weaker flows 
are partially truncated by the strong zonal flows of differential 
rotation. In the polar regions the convection is more isotropic 
and cyclonic. There the networks of downflow lanes surround 
upflows and both propagate in a retrograde fashion. 

The convection establishes a prominent differential rota- 
tion profile by redistributing angular momentum and entropy, 
building gradients in latitude of angular velocity and tempera- 
ture. Figures[TJ>,/ show the mean angular velocity £l(r,9) for 
cases D3 and H3, revealing a solar-like structure with a pro- 
grade (fast) equator and retrograde (slow) pole. Figures Ht,g 
present in turn radial cuts of 57 at selected latitudes, which 
are useful as we consider the angular velocity patterns real- 
ized here with faster rotation. These il(r, 9) profiles are aver- 
aged in azimuth (longitude) and time over a period of roughly 
200 days. Contours of constant angular velocity are aligned 
nearly on cylinders, influenced by the Taylor-Proudman theo- 
rem. 

In the Sun, helioseismology has revealed that the contours 
of angular velocity are aligned almost on radial lines rather 
than on cylinders. The tilt of 17 contours in the Sun may be due 
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FIG. 3. — Field line tracings of magnetic wreaths in case D3. (a) Snapshot of two wreaths in full volume at same instant as in Fig. [5] Lines trace the magnetic 
fields, color denoting the amplitude and polarity of the longitudinal field (red, positive; blue, negative). Magnetic field threads in and out of the wreaths, 
connecting the two opposite polarity structures across the equator (i.e., region A) and to the polar regions where the magnetic field is wound up by the cyclonic 
convection, (b) Same snapshot showing south polar region, (c) Zoom in on region A showing the complex interconnections across the equator between the two 
wreaths and to high latitudes. Convective flows create the distinctive waviness visible in all three images. 

to weakly couple the polar regions to the equatorial convec- 
tion. Flows of meridional circulation are slightly stronger 
in the dynamo cases than in the purely hydrodynamic cases, 
though both cases have weaker flows than are foun d in simula- 
tions r otating at the solar rate. Thus, as found in lBrown et al.l 
(2008), the flows of meridional circulation appear to weaken 
with more rapid rotation. The multi-celled nature of these 
meridional circulations may hold implication s for flux trans- 
port dynamo models ("e.g. lJouve & Brunl2007l) . Recent mean- 
field dynamo models are also beginning to explore the impli- 
cations of weaker and multi-celled meridional circ ulations for 
dynam o action in more rapidly rotating suns (e.g.. lJouve et al.1 
2009). 



TABLE 2 
Near-surface Af! 



Case 


AQia, 




Afiiat/fieq 


Epoch 


D3 


1.18 


0.71 


0.137 


2010-6980 


H3 


2.22 


0.94 


0.246 





NOTE. — Angular velocity shear in units of 
/xrads -1 , with Af2i at measured near the surface 
(0.97Rq) and AQ, measured across the full shell at 
the equator. The relative latitudinal shear Af\, t /f2 C q 
is also measured at the same point near the surface. 
For the dynamo case, these measurements are taken 
over the indicated range of days. Case D3 shows slow 
variations in Af2i at over periods of about 2000 days. 
The hydrodynamic case is averaged for roughly 300 
days and shows no systematic variation on longer 
timescales. 



quantified by measureme nts of the lat i tudina l angular velocity 
shear Afiiat. Here, as in lBrown et al.l ([2008), we define AQ\- dt 
as the shear near the surface between the equator and a high 
latitude, say ±60° 



Afilat = Hen — fifi 



(12) 



and the radial shear Ail r as the angular velocity shear between 
the surface and bottom of the convection zone near the equator 



97 R: 



-a 



0.72ft o • 



(13) 



We further define the relative shear as A51i at /f2 eq . In both 
definitions, we average the measurements of Afl in the north- 
ern and southern hemispheres, as the rotation profile is often 
slightly asymmetric about the equator. Case H3 achieves an 
absolute contrast A57i at of 2.22 \i rads -1 (352 nHz) and a rela- 
tive contrast of 0.247. The strong global-scale magnetic fields 
realized in the dynamo case D3 serve to diminish the differen- 
tial rotation. As such, this case achieves an absolute contrast 
Ar2i at of only 1.18 n rads" 1 (188 nHz) and a relative contrast 
of 0.137. This results from both a slowing of the equatorial 
rotation rate and an increase in the rotation rate in the polar 
regions. These results are quoted in Table|2] 

The meridional circulations realized in the dynamo case D3 
are very similar to those found in its hydrodynamic progeni- 
tor (case H3). As illustrated in Figures ld,h, the circulations 
are multi-celled in radius and latitude. The cells are strongly 
aligned with the rotation axis, though some flows along the in- 
ner and outer boundaries cross the tangent cylinder and serve 



3.2. Kinetic and Magnetic Energies 

Convection in these rapidly rotating dynamos is responsible 
for building the differential rotation and the magnetic fields. 
In a volume averaged sense, the energy contained in the mag- 
netic fields in case D3 is about 10% of the kinetic energy. 
About 35% of this kinetic energy is contained in the fluc- 
tuating convection (CKE) and about 65% in the differential 
rotation (DRKE), whereas the weaker meridional circulations 
contain only a small portion (MCKE). The magnetic energy is 
split between the contributions from fluctuating fields (FME), 
involving roughly 53% of the total magnetic energy, and the 
energy of the mean toroidal fields (TME) that are 43% of the 
total. The energy contained in the mean poloidal fields (PME) 
is only 4% of the total magnetic energy. These energies are 
defined as 



CKE=-p 



(v r - (vv)) 2 + (ve- (ve)) + 



■MY 



DRKE=-p(v ) 2 , 



MCKE 



(14) 

(15) 

(16) 
(17) 
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TABLE 3 
Energies 



Case 


CKE 


DRKE 


MCKE 


FME 


TME 


PME 


D3 


2.31 


4.35 


0.010 


0.36 


0.29 


0.029 


H3 


2.56 


22.2 


0.012 









NOTE. — Volume-averaged energy densities relative to the 
rotating coordinate system. Kinetic energies are shown for 
convection (CKE), differential rotation (DRKE) and merid- 
ional circulations (MCKE). Magnetic energies are shown for 
fluctuating magnetic fields (FME), mean toroidal fields (TME) 
and mean poloidal fields (PME). All energy densities are re- 
ported in units of 10 6 ergcirT 3 and are averaged over 1000 day 
periods. 



FME= ■ 



1 r 



(B r -{B r }) 2 +(Bg-(Bg)) 2 + 



TME= — (B^) 2 , 

OTT 

PME=-i- ((B r ) 2 +(B S ) : 



(18) 
(19) 
(20) 



where angle brackets denote an average in longitude. 

These results are in contrast to our previous simulations of 
the solar dynamo, where the mean fields contained only about 
2% of the m agnetic energy an d the fluctuating fields contained 
nearly 98% (Brun et al. 2004). In simulations of the solar dy- 
namo that inc luded a stable tachocl ine at the base of the con- 
vection zone dBrowning et al.ll2006l) . the energy of the mean 
fields in the tachocline can exceed the energy of the fluctuat- 
ing fields there by about a factor of three, though the fluctuat- 
ing fields still dominate the magnetic energy budget within the 
convection zone itself. Simulations of dynamo act ivity in the 
convecting cores of A-type stars dBrun et al.l 120051) achieved 
similar results. There in the stable radiative zone the energies 
of the mean fields were able to exceed the energy contained 
in the fluctuating fields, but in the convecting core the fluc- 
tuating fields contained roughly 95% of the magnetic energy. 
Simulations of dynamo action in fully-convective M-stars do 
howev er show high lev els of magnetic energy in the mean 
fields (Browning 2008). In those simulations the fluctuating 
fields still contain much of the magnetic energy, but the mean 
toroidal fields possess about 18% of the total throughout most 
of the stellar interior. In our rapidly rotating suns, the mean 
fields comprise a significant portion of the magnetic energy 
in the convection zone and are as important as the fluctuating 
fields. 

Convection is similarly strong in both rapidly rotating 
cases, and CKE is similar in magnitude. The differential rota- 
tion in the dynamo case is much weaker than in the hydrody- 
namic progenitor, and DRKE has decreased by about a factor 
of five. Meridional circulations are comparably weak in both 
cases. 

4. WREATHS OF MAGNETISM 

These rapidly rotating dynamos produce striking magnetic 
structures in the midst of their turbulent convection zones. 
The magnetic field is organized into large banded, wreath- 
like structures positioned near the equator and spanning the 
depth of the convection zone. These wreaths are shown for 



case D3 at two depths in the convection zone in Figure[2] The 
dominant component of the magnetic wreaths is the strong 
longitudinal field B$, with each wreath possessing its own po- 
larity. The average strength of the longitudinal field at mid- 
convection zone is ±7 kG and peak field strengths there reach 
roughly ±26 kG. Threaded throughout the wreaths are weaker 
radial and latitudinal magnetic fields, which connect the two 
structures across the equator and also to the high-latitude re- 
gions. 

These wreaths of magnetism survive despite being embed- 
ded in vigorous convective upflows and downflows. The 
convective flows leave their imprint on the magnetic struc- 
tures, with individual downflow lanes entraining the magnetic 
field, advecting it away, and stretching it into B, while leav- 
ing regions of locally reduced B^. The slower upflows carry 
stronger B^, up from the depths. Where the magnetism is 
particularly strong the convective flows are disrupted. Mean- 
while, where the convective flows are strongest, the longitu- 
dinal magnetic field is weakened and appears to vanish. In 
reality, the magnetic wreaths here are diving deeper below the 
mid-convection zone, apparently pumped down by the pum- 
meling action of the strong downflows. 

The deep structure of these wreaths is revealed by field line 
tracings throughout the volume, shown in Figure [3] for the 
same instant in time. The wreaths are topologically leaky 
structures, with magnetic field lines threading in and out of 
the surrounding convection. The wreaths are connected to the 
high-latitude (polar) convection, and on the poleward edges 
they show substantial winding from the highly vortical con- 
vection found there. This occurs in both the northern and 
southern hemispheres, as shown in two views at the same in- 
stant (north, Fig. [3jir and south, Fig. [3j>). It is here that the 
global-scale poloidal field is being regenerated by the cou- 
pling of fluctuating velocities and fluctuating fields. Mag- 
netic fields cross the equator, tying the two wreaths together 
at many locations (Fig. |3j:). The strongest convective down- 
flows leave their imprint on the wreaths as regions where the 
field lines are dragged down deeper into the convection zone, 
yielding a wavy appearance to the wreaths as a whole. 

4.1. Wreaths Persist for Long Epochs 

The wreaths of magnetism built in case D3 persist for long 
periods of time, with little change in strength and no rever- 
sals in global-scale polarity for as long as we have pursued 
these calculations. The long-term stability of the wreaths re- 
alized by the dynamo of case D3 is shown in Figure [4] Here 
the azimuthally-averaged longitudinal field (B^) and colatitu- 
dinal field (Bg) are shown at mid-convection zone at a point 
after the dynamo has equilibrated and for a period of roughly 
5000 days (i.e., several ohmic diffusion times). During this in- 
terval there is little change in either the amplitude or structure 
of the mean fields. This is despite the short overturn times of 
the convection (10-30 days) or the rotation period of the star 
(~ 9 days). The ohmic diffusion time at mid-convection zone 
is approximately 1300 days. 

Though the mean (global-scale) fields are roughly steady in 
nature (Figs. SJi, Z?), the magnetic field interacts strongly with 
the convection on smaller scales. Several samples of longitu- 
dinal field B^ are shown in full Mollweide projection at mid- 
convection zone (Fig. HJ:). The magnetic fields are clearly 
reacting on short time scales to the convection but the wreaths 
maintain their coherence. 

There are also some small but repeated variations in the 
global-scale magnetic fields. Visible in Figure [4j? are events 
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FIG. 4. — Persistent wreaths of magnetism in case D3. (a) Time-latitude plots of azimuthally-averaged longitudinal field (B^) at mid-convection zone (0.85/f v .; ) 
in a view spanning latitudes from ±70°, with scaling values indicated. The two wreaths of opposite polarity persist for more than 4000 days, (b) Mean 
colatitudinal magnetic field (Bg) at mid-convection zone over same interval, (c) Snapshots of B^ in Mollweide projection at mid-convection zone, shown for 
three times indicated in a,b. The wreaths maintain constant polarity over long time intervals, but still show variation as they interact with the convection. Time 
t2 corresponds to the snapshot in Fig.fSJ). 



where propagating structures of (Bg) reach toward higher lat- 
itudes over periods of about 1000 days (i.e., from day 3700 
to day 4500 and from day 5600 to day 6400). These are ac- 
companied by slight variations in the volume-averaged mag- 
netic energy densities and the comparable kinetic energy of 
the differential rotation. These variations are also visible in 
the differential rotation itself, as shown in Figure The dif- 
ferential rotation is fairly stable, though some time variation 
is visible at high latitudes. This is better revealed (Fig. |5p) 
by subtracting the time-averaged profile of il at each latitude, 
revealing the temporal variations about this mean. In the po- 
lar regions above ±40° latitude, speedup features move pole- 
ward over 500 day periods. These features track similar struc- 
tures visible in the mean magnetic fields (Fig.[4j?). The bands 
of velocity speedup bear some resemblance to the poleward 
branch of torsional oscillations observed in the solar convec- 
tion zone over the course o f a solar mag netic activity cycle 
(e.g., iThompson etaT]|2003t iHowd 120091 though here they 
propagate to higher latitudes on a shorter time scale. 

The temporal variations of the angular velocity contrast in 
latitude Af2i at are shown for this period in Figure [5J;. At mid- 
convection zone (sampled by red line) the variations in Af2i at 
are modest, varying by roughly 8%. Near the surface (green 
line) Af2i at shows similar variations with amplitudes of about 
6%. The near-surface values of Ai7i at are reported in Table|2] 
averaged over this entire period. 

These evolving structures of magnetism and faster and 
slower differential rotation appear to be the first indications 
of behavior where the mean fields themselves begin to wax 
and wane substantially in strength. As the magnetic Reynolds 
number is increased, by either decreasing the magnetic diffu- 
sivity r\ or by increasing the rotation rate of the star fio, this 
time varying behavior becomes more prominent and can even 
result in organized changes in the global-scale polarity. Such 
behavior is evident in a number of our dynamo simulations 
and will be reported on in a subsequent paper. 

5. CREATING MAGNETIC WREATHS 



tie variation in time. To understand the physical processes 
responsible for maintaining these magnetic wreaths, we ex- 
amine the terms arising in the time- and azimuth-averaged in- 
duction equation for case D3. 

5.1. Maintaining Wreaths of Toroidal Field 

We begin our analysis by exploring the maintenance of the 
mean toroidal field (B$). Here it is helpful to break the in- 
duction term from equation into contributions from shear, 
advection and compression, namely 



V X (v X B) = 

(£^V)w-(^y)S- B(V-v) . 

shear advection compression 



(21) 



Details of this decomposition are given in the Appendix. 

The evolution of the mean longitudin al (to roidal) field (B^) 
is described symbolically in e quation dA8b . with individual 
terms defined in equation dA9l ). When we analyze these terms 
in case D3, we find that (Bj,) is produced by the shear of dif- 
ferential rotation and is dissipated by a combination of tur- 
bulent induction and ohmic diffusion. This balance can be 
restated as 



dt 



' Pvis + (Pfs + Pfa + Pmd) « , (22) 



with Pms representing production by the mean shearing flow 
of differential rotation, Pps by fluctuating shear, Pfa by fluc- 
tuating advection, and Pmd by mean ohmic diffusion. Those 
terms are in turn 



Pms= ((B)-V)(v)\^ 
P FS = <(5'.V),%, 
Pfa = -(K-V)5%, 

Pmd = -V X77V X (B)\ 



(23) 
(24) 
(25) 
(26) 



The magnetic wreaths formed in case D3 are dominated 
by strong mean longitudinal field components and show lit- 



where brackets again indicate an azimuthal average and 
primes indicate fluctuating terms: v' = V — (v). The detailed 
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FIG. 5. — Differential rotation in case D3. (a) Angular velocity Q at mid- 
convection zone (0.85Sq), with ranges in both nHz and /xrads~' . The equator 
is fast while the poles rotate more slowly, (b) Temporal variations are empha- 
sized by subtracting the time-averaged profile of f2(r, ff), revealing speedup 
structures at high latitudes and pulses of f ast a nd slow motion near the equa- 
tor, (c) Angular velocity shear Af2i at (eg. 1121 near the surface (upper curve, 
green) and at mid-convection zone (lower, red). 

implementation of the se terms is presented for our spherical 
geometry in equations ( IA10HAl5T l. These terms are illustrated 
in Figure[6]for case D3, averaged over a 450 day interval from 
day 6450 to 6900. 

The structure of (B^) is shown in Figure^. The shearing 
flows of differential rotation P M s (Fig. |6£>) act almost every- 
where to reinforce the mean toroidal field. Thus the polarity of 
this production term generally matches that of (B^). This pro- 
duction is balanced by destruction of mean field arising from 
both turbulent induction and ohmic diffusion (sum shown in 
Fig-ISt). The individual profiles of Pps, ^fa and Pmd are pre- 
sented in turn in Figures |6j/,e,/. The terms from turbulent 
induction (Pps and Pfa) contribute to roughly half of the to- 
tal balance, with the remainder carried by ohmic diffusion of 
the mean fields (Pmd)- In the core of the wreaths, removal of 
mean toroidal field is largely accomplished by fluctuating ad- 
vectionPpA (Fig.|6^) and mean ohmic diffusion Pmd (Fig.|6]0, 
with the latter also important near the upper boundary. Turbu- 
lent shear becomes strongest near the bottom of the convec- 
tion zone and in the regions near the high-latitude side of each 
wreath. Thus Pps (Fig.|6]i) becomes the dominant member of 
the triad of terms seeking to diminish the mean toroidal field 
there. We find that the mean poloidal field is regenerated in 



roughly the same region. 

In the analysis presented in Figure|6]we have neglected the 
advection of (B$) by the meridional circulations (shown in 
the Appendix as Pma), which we find plays a very small role 
in the overall balance. We have also neglected the amplifi- 
cation of (Btp) by compressibility effects (the Appendix, Pmc 
and PfcX though it does contribute slightly to reinforcing the 
underlying mean fields within the wreaths. 

To summarize, the mean toroidal fields are built through 
an J7-effect, where production by the mean shearing flow of 
differential rotation (Pus) builds the underlying (B^). In the 
statistically steady state achieved, this production is balanced 
by a combination of turbulent induction (Pps +Pfa) and ohmic 
diffusion of the mean fields (Pmd)- 

5.2. Maintaining the Poloidal Field 

The production of mean poloidal field is achieved through a 
slightly different balance, with turbulent induction producing 
poloidal field and ohmic diffusion acting to dissipate it. The 
mean flows play little role in the overall balance. This balance 
is clarified if we represent the mean poloidal field by its vector 
potential (A^), where 



(27) 



(B pol ) = (B l .)r+(B e )e = Vx{A 4> cf> 



as discussed in the Appendix. We recast the induction equa- 
tion in terms of the poloidal vector potential by uncurling 
the equation once, obtaining 



dt 



{vxB}\,- v Vx(B)\ 



(28) 



which is also equation dA29l ) in the Appendix. The first term 
is the electromotive force (emf) arising from the coupling of 
flows and magnetic fields, and the second term is the ohmic 
diffusion. These can be decomposed into contributions from 
mean an d fluct uating components, as shown symbolically in 
equation dA30b . 

In case D3 we find that the mean poloidal vector potential 
(A^) is produced by the fluctuating (turbulent) emf and is dis- 
sipated by ohmic diffusion 



dt 



! Ern + E\ 



0. 



(29) 



with £pi the emf arising from fluctuating flows and fluctuating 
fields, and contributing to the mean induction. The Emd is the 
emf arising from mean ohmic diffusion. These terms are 



E P1 =(v'xB% = (v' r B> ei 
£md = -?7Vx (B)\,. 



(30) 
(31) 



The co ntribu tion arising from the omitted term Eya 
(see eq. IA31I ). related to the emf of mean flows and mean 
fields, is smaller than these first two by more than an order of 
magnitude. Additionally, Eya has a complicated spatial struc- 
ture which does not appear to act in a coherent fashion within 
the wreaths to either build or destroy mean poloidal field. 

The mean vector potential (A^) is shown in Figure|7Ji, with 
poloidal field lines represented by the overlying contours. The 
mean radial magnetic field (B r ) is about ±1 kG in the cores of 
the wreaths, whereas the mean colatitudinal field (Bg) has an 
amplitude of roughly -2 kG (thus directed northward in both 
hemispheres), concentrated near the bottom of the convection 
zone. 
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FIG. 6. — Generation of mean toroidal magnetic field in case D3. The view is from ±45° latitude to emphasize the equatorial regions, (a) Mean toroidal field 
(B^) with wreaths strongly evident, (b) Production by Pms serves to build (B$). This rate term generally matches the sense of (B^), thus being negative (blue 
in colorbar, with ranges indicated) in the core of the northern wreath and positive (red) in that of the southern wreath, (c) Destruction of mean toroidal field is 
achieved by the sum of the two fluctuating (turbulent) induction terms and the ohmic diffusion (Pps +Pfa +Pmd)- This sum clearly has opposite sense and similar 
magnitude to Pms- We break out these three destruction terms in the following panels, (d) Fluctuating (turbulent) shear Pps is strongest near the high-latitude 
side of each wreath, and (e) fluctuating (turbulent) advection Pfa is strongest in the cores of the wreaths. The sum of these terms (Pps +Pfa) is responsible for 
about half the destructive balance, with the remainder coming from (/) the mean ohmic diffusion Pmd- Some differences arise in the boundary layers at top and 
bottom. 

field is here being converted into mean poloidal field by the 
fluctuating flows. 

There are two terms that contribute to Epi, as shown in 
equation (f30b . Much of that fluctuating emf arises from corre- 
lations between fluctuating latitudinal flows and radial fields 
(—v' e B' r ), which follows the structure of Epi (Fig.[7J?) closely. 
The contribution from fluctuating radial flows and colatitudi- 
nal fields (v' r B' s ) is more complex in structure. Near ±20° 
latitude, this term reinforces (— v' e B' r ), but acts against it at 
higher latitudes and thus diminishes the overall amplitude of 
£fi- The mean ohmic diffusion £md (Fig.|7K), almost entirely 
balances the production of (A^) by Epi. 

This shows that our mean poloidal magnetic field is main- 
tained by the fluctuating (turbulent) emf and is destroyed by 
ohmic diffusion. In mean-field dynamo theory, this is often 
parametrized by an "a-effect." Now we turn to interpretations 
within that framework. 

6. EXPLORING MEAN-FIELD INTERPRETATIONS 

Many mean-field theories assert that the production of mean 
poloidal field is likely to arise from the fluctuating emf. This 
process is often approximated with an a-effect, where it is 
proposed that the sense and amplitude of the emf scales with 
the mean toroidal field 




FIG. 7. — Production of mean poloidal vector potential (As) in case D3, 
with view restricted to ±45° latitude, (a) Mean poloidal vector potential 
(A*), with sense denoted by color (red, clockwise; blue, counter-clockwise). 
(b) The fluctuating (turbulent) emf £pi acts to build the vector potential. This 
term is strongest near the bottom of the convection zone and the poleward side 
of the wreaths, (c) Mean ohmic diffusion £"md acts everywhere in opposition 
to Em. The cores of the wreaths are positioned at roughly ±15° latitude 
(Fig.ll/). 

The production of (Aj,) by the fluctuating (turbulent) emf 
Epi is shown in Figure[7I?. Here too we average over the same 
450 day interval. This term generally acts to reinforce the 
existing poloidal field, having the same sense as the under- 
lying vector potential in most regions. It is strongest near the 
bottom of the convection zone and is concentrated at the pole- 
ward side of each wreath. This is similar, though not identi- 
cal, to the structure of destruction of mean toroidal field by 
fluctuating shear P F s (Fig. [6}/). It suggests that mean toroidal 



(v'xB')=a(B) 



(32) 



where a can be either a simple scalar or may be related to 
the kinetic and magnetic (current) helicities. In isotropic (but 
not reflectionally symmetric), homogeneous, incompressible 
MHD turbulence 



a= -(a k + a m ), 
ak = -v' ■ (V X v') , 
a m = ^B'-(VXB'), 



(33) 
(34) 
(35) 
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FIG. 8. — Estimating the mean-field ct-effect from case D3. Shown are 
the ( a) kine tic and (b) magnetic contributions to the a-effect as defined in 
eqs. 1331351 . (c) Mean-field a, constructed by combining and q,„ with a 
turbulent correlation time r. 

as discussed in iPouauet et al] (119761) and 
iBrandenburg & Subramanianl (120051) . Here r is the life- 
time or correlation time of a typical turbulent eddy. In 
mean-field theory, these fluctuating helicities are typically 
not solved directly and are instead solved through auxiliary 
equations for the total magnetic helicity or are prescribed. 
Here we can directly measure our fluctuating helicities and 
examine whether they approximate our fluctuating emf. 

To assess the possible role of an a-effect in our simulation, 
we show in Figures \&fr,b the fluctuating kinetic and current 
helicities a* and a m realized in our case D3, averaged over 
the same 450 day analysis interval. To make an estimate of the 
a-effect, we approximate the correlation time r by defining 



where H P is the local pressure scale height and v' is the local 
fluctuating rms velocity, which are functions of radius only. 
Estimated by this method, the turnover time t has a smooth 
radial profile and is roughly 10 days near the bottom of the 
convection zone, 3 days at mid-convection zone, and slightly 
less near the upper boundary. If we use the fast peak upflow 
or downflow velocities instead of the rms velocities, our es- 
timate of r is about a factor of 4 smaller. Our mean-field a 
(eq.l33Tl is shown in Figure[8];. In the upper convection zone, 
this is dominated by the fluctuating kinetic helicity while the 
fluctuating magnetic (current) helicity becomes important at 
depth. 

We form a mean-field emf (right-hand side of eq. [32| > by 
multiplying our derived a (Fig. [§}:) with our (B^,) (Fig. |6ji), 
and show this in Figure |9ji. The turbulent emf E-pi, which is 
the left-hand side of equation (l32l l. can be measured in our 
simulations and is shown again in Figure [9j?. Although there 
is some correspondence in the two patterns, there are signifi- 
cant differences. In particular, the mean-field emf a(B ( p) has 
peak amplitudes in the cores of the wreaths (at ±15° latitude) 
and is negative there. In contrast, the actual fluctuating emf 
given by £pi is positive and has its highest amplitude at the 
poleward side of the wreaths (near ±20° latitude). Thus the 
mean-field emf predicts an incorrect balance in the genera- 




FlG. 9. — Comparison of emfs in case D3. (a) Profile of proposed mean- 
field emf given by ct(Bj,). (b) Actual turbulent emf £pi measured in the 
dynamo, (c) Variation of hemisphere-averaged emfs with fractional radius. 
The mean-field approximated emf is shown in blue, and £pi in red. The 
average over the northern hemisphere is shown solid, the southern is dashed. 



tion terms and would yield a distinctly different mean poloidal 
magnetic field. 

To assess whether better agreement may be achieved with 
a latitude-averaged emf, we average the mean-field emf and 
Epi separately over the northern and southern hemispheres 
and plot these quantities in Figure [9}:. Though both have a 
similar positive sense near the base of the convection zone, 
the hemisphere-averaged £fi becomes small above O.SRq 
whereas the averaged mean-field emf a(B ( j > ) is large and neg- 
ative there. Thus even the averaged emfs are not in accord. 

In summary, it is evident that a simple scalar a-effect will 
predict the wrong sign for the fluctuating emf in the two 
hemispheres, as (B^) is anti-symmetric across the equator 
while (A^) is symmetric. An a-effect based on the kinetic 
helicity and magnetic helicity may capture some sense of 
the fluctuating emf, as those quantities are themselves anti- 
symmetric across the equator. Yet Figure|9]suggests that there 
are significant discrepancies between this particular approxi- 
mation and our turbulent emf. In particular, this mean-field 
a-effect misses the offset between the generation regions for 
mean toroidal and mean poloidal field. This offset in lati- 
tude of the generation regions may be important for avoid- 
ing the a-quenching problems encountered in many mean- 
field theories. A more complex mean-field model, which 
takes spatial gradients of (B$) into acc ount, may do bet- 
ter. I n particular, the fl X ./-effect (e .g., iMoffatt & Proctor! 
119821: Rogachevski i & Kleeorinl 120031) may be at work in 
these systems, and preliminary explorations indicate that this 
term matches the spatial structure of our £fi better than the 
above a-effect. A tensor representation of the a-effect may 
also do much better at approximating E Fl , and test-field tech- 
niques could be em ployed to measure this quantity (e.g., 
Schrinner et al. 2005, and recently reviewed in Brandenburg 
120091) . As with our analysis of dynamo production terms pre- 
sented in $5] this comparative study of a(B ( p) and Epi is con- 
ducted here for the special circumstances of a dynamo which 
builds global-scale magnetic fields that are nearly steady in 
time. The magnetic wreaths realized in dynamos at higher 
magnetic Reynolds numbers show larger time-variations, and 
it is possible that a(B t p) better approximates £pi during the 
growing phase of each oscillation, when the magnetic fields 
have not yet saturated in strength and the dynamo is in a more 
kinematic regime. 
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7. CONCLUSIONS 

The ability for a dynamo to build wreaths of strong mag- 
netic fields in the bulk of the convection zone has largely 
been a surprise, for it had generally been supposed that tur- 
bulent convection would disrupt such magnetic structures. 
To avoid these difficulties, many solar and stellar dynamo 
theories shift the burden of magnetic storage, amplification 
and organization to a tachocline of shear and penetration at 
the base of the convection zone where motions are more 
quiescent. In contrast, our simulations of rapidly rotating 
stars are able to achieve sustained global-scale dynamo action 
within the convection zone itself, with the magnetic struc- 
tures both being built and able to survive while embedded 
deep within the turbulence. These dynamos are able to cir- 
cumvent the Parker instability by means of turbulent Reynolds 
and Maxwell stresses that contribute to the mechanical force 
balance and prevent the wreaths from buoyantly escaping the 
convection zone. This striking behavior may be enabled by 
the stars rotating somewhat faster than the current Sun, which 
yields a strong differential rotation that is a key element in the 
dynamo behavior. In our broader exploration of rapidly ro- 
tating dynamos, we find that magnetic wreaths are present in 
all simulations, including those rotating as slowly as 1.5^0. 
Such structures may be obtainable in simulations rotating at 
the solar rate as well, and efforts are underway to explore the 
presence of wreaths in solar dynamos. 

We have achieved some dynamo states that are persistent 
and others that flip the sense of their magnetic fields. In our 
case D3 the global-scale fields have small vacillations in their 
amplitudes, but the magnetic wreaths retain their identities for 
many thousands of days. This represents hundreds of rotation 
periods and several magnetic diffusion times, indicating that 
the dynamo has achieved a persistent equilibrium. 

Increasing the rotation rate or decreasing the magnetic dif- 
fusivity rj yields more complex time dependence. In many 
of our dynamos the oscillations can become large, and this 
may result in the global-scale fields repeatedly flipping their 
polarity. At times those dynamos appear to be cyclic but 
in other intervals they behave more chaotically. Such time- 
dependent dynamos will be reported on in a forthcoming pa- 
per. In separate explorations, we have found that magnetic 
wreaths also survive in the presence of a tachocline of pene- 
tration and shear. In those simulations the wreaths continue 
to fill the convection zone even while developing roots in the 
tachocline. Dynamos in rapidly rotating suns with tachoclines 
can also exhibit time-dependent oscillations and polarity re- 
versals. Wreath-building dynamos with tachoclines will be 
reported on subsequently. 

In our persistent case D3 we are able to analyze the gener- 
ation and transport of mean magnetic field. We find that our 
dynamo action is of an a - £1 nature, with the mean toroidal 
fields being generated by an £!-effect from the mean shear- 
ing flow of differential rotation. This generation is balanced 
by a combination of turbulent induction and ohmic diffusion. 
The mean poloidal fields appear to be generated by an a- 
effect arising from couplings between the fluctuating flows 
and fluctuating fields, with this production largely balanced 



by the ohmic diffusion. This is unlike the toroidal balance, 
for here the mean flows play almost no role and the turbu- 
lent correlations are constructive rather than destructive. In 
assessing what a mean-field model might predict for the mag- 
netic structures realized in case D3, we find that the isotropic, 
homogeneous a-effect based on kinetic and magnetic (cur- 
rent) helicities fails to capture the sense of our turbulent emf. 
In general, our £pi is poorly represented by an a{B$) that is 
so determined. This comparative analysis of a{B^) and £pi 
is performed here only for the special case of a dynamo with 
persistent global-scale magnetic fields. It is possible that these 
results will differ in our dynamos that show substantial time- 
varying behavior. 

The realization of global-scale magnetic structures in our 
simulations, and their great strength relative to the fluctuating 
fields, may in part be a consequence of the relatively modest 
degree of turbulence attained here. Whether such structures 
can be generated and sustained amidst the far more complex 
flows in actual stellar interiors is not yet clear. If such struc- 
tures are indeed realized in stars, they may or may not sur- 
vive to print through the highly turbulent convection occur- 
ring just below the stellar photosphere. If they do appear at 
the surface, some global-scale magnetic features may propa- 
gate toward the poles along with the bands of angular velocity 
speedup. There are some indications in stellar observations 
that global-scale toroidal magn etic fields may indeed become 
strong in rapidly rotating stars dDonati et al.ll2006t iPetit et all 
2008), though small-scale fields may still account for much o f 
the magnetic energy near the surface dReiners & Basrill2009l) . 
The global-scale poloidal fields may be more successful in 
surviving the passage through the turbulent surface convec- 
tion. If they do, the stellar magnetic field will likely have 
significant non-dipole components. Thus the mean poloidal 
fields observed at the surface may give clues to the presence 
of large wreaths of magnetism that occupy the bulk of the 
convection zone. 
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APPENDIX 

PRODUCTION, DESTRUCTION AND TRANSPORT OF MAGNETIC FIELD 

We derive diagnostic tools to evaluate the generation and transport of magnetic field in a magnetized and rotating turbulent 
convection zone. This derivation is in spherical coordinates, and is under the anelastic approximation. 
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Induction Equation 

In the induction equation ©, the first term on the right hand side represents production of magnetic field while the second 
term represents its diffusion. We first rewrite the production term to make the contributions of shear, advection and compressible 
effects more explicit as 

V x(vxB) = (B- V)v-(v ■ V)B-B(V-v). (Al) 

Under the anelastic approximation the divergence of v can be expressed in terms of the logarithmic derivative of the mean density 
because 

V- (pv) = = p( V- v) + (v ■ V)p, 

and therefore 

„ d 

V-v=-v r —lnp. (A2) 
or 

The induction equation thus becomes 

— = (B-V)v-(v V)B + v r B—\np - V X (r]V X B) (A3) 

shearing advection v * . diffusion 
& compression 

As labeled, the first term represents shearing of B, the second term advection of B, the third one compressible amplification of 
B, and the last term ohmic diffusion. 



Production of Axisymmetric Magnetic Field 

To identify the processes contributing to the production of mean (axisymmetric) field, we separate our velocities and magnetic 
fields into mean and fluctuating components v = (v) + v' and B = (B) + B' where angle brackets denote an average in longitude. 
Thus (v') = (B 1 ) = by definition. Expanding the production term of equation (I A3t we obtain the mean shearing term 



the mean advection term 



and the mean compressibility term 



((B ■ V)v) = ((B) ■ V) («> + ({B' ■ V)u' >, (A4) 
- ((« • V)B) =-((v)- V) (B) - ({v 1 ■ V)B'), (A5) 



(v r B-\np) = ((v r ){B) + (v' r B'}) -hip. (A6) 
or v 'Or 



In a similar fashion, the mean diffusion term becomes 

-(Vx(ijVxB))=-Vx(r)Vx(B)). (A7) 
The axisymmetric component of the induction equation is written symbolically as: 

1 ~ = ^ms + Pfs + Pma + Pfa + Pmc + Pfc + Pmd ( A8) 

ot 

With Pms representing production of field by mean shear, Pps production by fluctuating shear, Pma advection by mean flows, fpA 
advection by fluctuating flows, Pmc amplification arising from the compressibility of mean flows, Ppc amplification arising from 
fluctuating compressible motions, and Pmd ohmic diffusion of the mean fields. In turn, these terms are 

Pms = ((B)-V)(v), P fs =((B'-V)v'), P ma =-((v)-V){B), P fa =-((v'-V)B'), 

d d 

Pmc = ((v,.)(B))— lnp, P PC = ({v' r B'))— Inland P MD =- V X (r/V X (B)). (A9) 
We now expand each of these terms into their full representation in spherical coordinates. 
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dt 

Pms = 



-PmS + PfS + PmA + ^FA + PmC + ^FC + ^MD 

{B 4 )(v r ) + coiO(B 4 ){vg) 



Pfs - 
Pma = ~ 
Pfa = ~ 



B 



,d_ + Kd_ A B '<t> g 



' dr r dO rsin 



(BW r ) + cot9(BW e ) 



K rl 8r r 89 



(B„ 



(v )(fl r ) + cot0(v )(B ) 



,8 v' d v'a 
v' r — + - l — + - 



8 



" dr r d9 rsin defy 
d 



B'a 



(v'B' r ) + cotO(v'B' e ) 



Puc= «v r )(B ))-lnp 



Pmd= r,V 2 (B, 



dr 

r){B 4 ) dr] ( 1 d(r(B 4 )) 



r 2 sin 2 6> dr \r dr 
Production of Mean Latitudinal Field 



(A10) 
(All) 
(A12) 

(A13) 
(A14) 
(A15) 



d(B ) 



dt 
Pms = 

Pfs = 

Pma = ~ 

Pfa = ~ 



= Pms + Pfs + Pma + Pfa + Pmc + Pfc + Pmd 



dr r d6 



{ve) + 



(Bg)(v r )-cotO(B 4> )(v 4 ) 



,d B' d B' d 



' dr r d9 r sin d(j> 



(fi;,.;)-cot9{«>;) 



dr r dO 

, d v'a d 
v' — + — — + ■ 



(Be)~ 
v's d 



{v e ){B r )-cotO{v 4> ){B 4> ) 



" dr r dO rsin 6* d(f> 



r 



B'e)- 



(v> e B>)-c tO(v>^} 



Puc= {(v r )(B 6 ))-\np Pp C = ((v' r B' e )) -Inp 



Pmd = vV 2 (B e } + 



2r)d(B r ) ri{B g } dr, ( \ d(r{B e )) 1 d(B r ) 



r 2 dO r 2 sin 2 dr \r dr r dO 
Production of Mean Radial Field 



(A16) 
(A17) 
(A18) 
(A19) 
(A20) 
(A21) 



8(B r ) 
dt 

Pms = 
Pfs = 
Pma = ~ 
Pfa=~ 



Pms + ^fs + ^ma + ^fa + Pmc + Pfc + Pmd 

. (Bg)(vg) + {B 4> ){V 4 



IB \ 9 + {Be) 9 
{Br) d~r + — dO 

B ' 9B ± 9 . B * 9 



' dr r dO rsin d<f> 



(B g v' e ) + (B'v') 



3 Mi 

W 5r r dO 
t 9 , v 'e d , 



(B r ) + 
v's d 



(vg)(Bg) + ( V<p ){B 4 



" dr r dO rsin# d<\> 



b:)+ 



Pmc= {{v r ){B r ))^-\np P FC =((vX))|-lnp 



MD : 



v V 2 (B r )-2 V ^- 



(B r ) 2r/ d(B g ) 277 cot (Bg) 



r 2 dO 



(A22) 
(A23) 
(A24) 
(A25) 
(A26) 
(A27) 
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Maintaining the Poloidal Vector Potential 

The balances achieved in maintaining the mean poloidal magnetic field are somewhat clearer if we consider its vector potential 
rather than the fields themselves. The mean poloidal field (B po \) has a corresponding vector potential (A^), where 

(B poi ) =(B r )r + (B e )0 = Vx(A\ 4> ) 

^M^^^-^K)^ ( A28) 

= v x (a^ 

The other components of the poloidal vector potential disappear, as terms involving d/d(f) vanish in the azimuthally-averaged 
equations. Likewise, the (^-component of the possible gauge term VA is zero by virtue of axisymmetry. We recast the induction 
equation (eq. |4j in terms of the poloidal vector potential by uncurling the equation once and obtain 

^M= vX B\^- V VxB\^. (A29) 
This can then be decomposed into mean and fluctuating contributions, and represented symbolically as 



dt 



= E M i + E F i + E M d, (A30) 



with Emi representing the electromotive forces (emf) arising from mean flows and mean fields, and related to their mean induction. 
Likewise, £pi is the emf from fluctuating flows and fields and Emd is the emf arising from mean diffusion. These are in turn 

E m = {v)x(B)\ <j> =(v r }(B e }-(vg}(B r ), (A31) 

E m = (v'xB')\^=(v' r B' g )-(v' e B' r ), (A32) 

E MD = - V V X {B)\,=-v\ {jM B e))-^§f) (A33) 

Fluctuating (Non-Axisymmetric) Component of the Induction Equation 

Left out of this analysis is the fluctuating component of the induction equation. This can be derived by subtracting the mean 
induction equation (IA8b from the full induction equation, yielding the following equation for the fluctuating fields 

BR' 

— = ({B)-V)v' + (B'-V){v)+e 
-((v)-V)B'-(v'-V)(B)-T 

+((v r )B' + v' r (B))^-lnp+g 
or 

-V X (77V X (B')) (A34) 

where the quantities £ = (B' ■ V)v' - ((B' ■ V)v'), T = (V • V)B'-((v' ■ V)B'), and Q = (v' r B' - {v' r B'))§-\np, represent the 
difference between mixed stresses from which we subtract their axisymmetric mean. In the standard mean-field derivation, these 
quantities are siblings of the G -current involv ing the mean electromotive force (v X B) and its 3-D equivalent v X B (i.e., the so 
called "pain in the neck" term, Motfanll 1978b. 
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